Boson-Fermion coherence in a spherically symmetric harmonic trap 
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We consider the photoassociation of a low-density gas of quantum-degenerate trapped fermionic 
atoms into bosonic molecules in a spherically symmetric harmonic potential. For a dilute system and 
the photoassociation coupling energy small compared to the level separation of the trap, only those 
fermions in the single shell with Fermi energy are coupled to the bosonic molecular field. Introducing 
a collective pseudo-spin operator formalism we show that this system can then be mapped onto the 
Tavis-Cummings Hamiltonian of quantum optics, with an additional pairing interaction. By exact 
diagonalization of the Hamiltonian, we examine the ground state and low excitations of the Bose- 
Fermi system, and study the dynamics of the coherent coupling between atoms and molecules. In a 
semiclassical description of the system, the pairing interaction between fermions is shown to result 
in a self-trapping transition in the photoassociation, with a sudden suppression of the coherent 
oscillations between atoms and molecules. We also show that the full quantum dynamics of the 
system is dominated by quantum fluctuations in the vicinity of the self-trapping solution. 

PACS numbers: 03.75.Lm, 03.75. Ss, 42.50.Ar 



I. INTRODUCTION 

The formation of ultracold diatomic molecules from 
Feshbach resonances and photoassociation has witnessed 
spectacular developments in recent years. Early demon- 
strations of molecule formation using two-photon Ra- 
man photoassociation and a Feshbach resonance 
were dominated by the molecular losses due to processes 
such as inelastic decay to lower energy molecular vibra- 
tional states Q, so that the existence of the molecules 
could only be inferred from the decrease in the num- 
ber of atoms. The first unambiguous coherent conversion 
of atoms into molecules was performed by Donley et al. 
P|, who exploited a Feshbach resonance in a 85 Rb Bose- 
Einstein Condensate (BEC). In subsequent experiments 
starting from an atomic condensate of 87 Rb, Rempe and 
coworkers used adiabatic rapid passage to create the 
molecules |||. Because the molecules and atoms have 
different magnetic moments, they could be spatially sep- 
arated from each other using a magnetic field gradient 
via the Stern-Gerlach effect. Similar work has been con- 
ducted by Xu et al. @- Starting from a sodium BEC, 
they used resonant laser light to blast away the remaining 
atoms in the sample and isolated the molecules. Unfortu- 
nately, the conversion efficiency was limited by inelastic 
losses very close to the resonance so that molecular yields 
were < 10%. 

For fermionic atoms close to a Feshbach resonance the 
inelastic collision rate for relaxation to lower energy vi- 
brational states of the molecules scales like a(B)~ 2 ' 55 
whereas for bosons it scales like a(B) where a(B) is the 
scattering length near the resonance Q ■ This is because 
close to resonance the effective size of the molecules is 
of the order a(B), which is comparable to the interpar- 
ticle spacing. In order for a molecule to decay to a more 
deeply bound vibrational state with radius R e <C a(B), 
the atoms comprising the molecule along with an ad- 
ditional atom must all collide within a distance ~ R e . 



Since two of the three atoms are necessarily identical, the 
collision rate is suppressed for fermions. Taking advan- 
tage of this consequence of the Pauli Exclusion Principle, 
Greiner et al. jij achieved the first molecular BEC by 
starting from a quantum degenerate spin mixture of 40 K 
using adiabatic rapid passage through a Feshbach reso- 
nance with a conversion efficiency of ~ 80%. At about 
the same time Zwierlein et al. and Jochim et al. [T(il ] 
succeeded in producing a BEC of 6 Li 2 dimers by evap- 
oratively cooling the atoms at a constant magnetic field 
just below a resonance where a(B) is large and positive. 
Two recent experiments have led to the observation of 
heteronuclear Feshbach resonances in Bose-Fermi mix- 
tures of 6 Li and 23 Na 0] m one case, and of 87 Rb and 
40 K in the other ^^|- We also mention recent experi- 
ments by Kerman et al. |l3[ . who produced metastable 
RbCs molecules in their lowest triplet state starting from 
a laser-cooled mixture of 85 Rb and 133 Cs by photoas- 
sociation. Major current experimental and experimen- 
tal efforts are directed towards the exploration of the 
crossover between the Bardeen-Cooper-Schrieffer (BCS) 
state of fermionic atoms and BEC of bosonic molecules 
as the system crosses a Feshbach resonance 0|. In the 
vicinity of these resonances, the system enters a strongly 
interacting regime that offers a challenge for many-body 
theories |l5l Il6|. 

It is known that the pairing properties of finite-size 
systems can be significantly different from those of the 
bulk material, due to the discrete energy spectrum of the 
particles involved. The detailed role of the shell struc- 
ture has been explored both in the nuclei 0] and of 
superconductor grains 18], in which a collective char- 
acter of pairs plays an important role. In this paper 
we consider the photoassociation of a dilute quantum- 
degenerate gas of fermionic atoms trapped in a spher- 
ically symmetric harmonic trap into molecular dimers, 
including the pairing interaction between fermions. By 
introducing a pseudo-spin formalism for the time rever- 
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sal pairing operator |19ll20| . this problem can be mapped 
onto an extension of the Tavis-Cummings model |2ll |22| 
that describes the coupling of N two-level atoms and a 
single mode of the electromagnetic field, and has found 
applications in the study of superradiance in quantum 
optics [23L l24j . This analogy allows us to study in detail 
the ground state and lower excited states of the system, 
as well as the coherent dynamics of atom-molecule cou- 
pling in the trap. We show that for appropriate condi- 
tions, only those fermions on the last energy shell of the 
trap participate in the photoassociation process, and dis- 
cuss the impact of the filling of that shell on molecule for- 
mation. We find that depending on the detuning S of the 
photoassociation laser from the energy difference between 
the molecules and atom pairs, the nature of the ground 
state changes from being predominantly atomic to pre- 
dominantly molecular in nature. We study the crossover 
between these two regions in detail, and quantify its prop- 
erty via the joint coherence of the atomic and molecular 
fields and the entanglement entropy of the system. 

The coherent conversion of fermions into bosons has 
been studied in the homogeneous case by several au- 
thors 0, |2^ |2|j. Trapped systems, in addition to 
being closer to the experimental situation, present sev- 
eral unique characteristics: First, the discreteness of the 
energy levels eliminates many of the difficulties associ- 
ated with a continuum. In addition, the high degener- 
acy of spherically symmetric harmonic potentials simpli- 
fies significantly the study of coherent quantum dynam- 
ics. Indeed, the problem resembles then the dynamics 
of a bosonic Josephson Junction j23, HI], although the 
nonlinear coupling between fermionic atoms and bosonic 
molecules leads to considerably richer dynamics. More- 
over the additional pairing interaction between fermions 
is shown to result in a self-trapping transition [27L l2^ | , 
with a sudden suppression of the coherent oscillations. 

This paper is organized as follows. Section II discusses 
our model and formulates it in terms of a pseudo-spin 
formalism. Section III presents results of the static prob- 
lem where the many-body states are classified by number 
of unpaired fermions known as seniority [sH ] in nuclear 
physics. We show that for a fixed number of atoms, the 
ground state of the system is always the state of mini- 
mum seniority. We examine the crossover behavior of the 
ground state as the detuning parameter 5 is varied, as a 
function of the ratio of the total number of fermionic pairs 
and molecules to the degenerate number of the Fermi 
level. The entanglement between fermions and bosons is 
evaluated and found to be reduced as the pairing inter- 
action becomes stronger. In section IV, we analyze the 
coherent dynamics of the nonlinear atom-molecule cou- 
pling. Using a semiclassical factorization ansatz, we show 
the appearance of a self-trapping transition in the pres- 
ence of pairing interaction. An exact quantum solution 
shows that around that transition point the dynamics 
is characterized by large quantum fluctuations. Finally, 
section V is a summary and outlook. Calculational de- 
tails are relegated to an appendix. 



II. MODEL 

We consider a trapped dilute gas of two-component 
fermionic atoms in hyperfine states of spin a ="f , J, at zero 
temperature and coupled to a single-mode gas of bosonic 
molecules via a two-photon Raman transition. The trap 
is assumed to be harmonic and spherically symmetric, 
described by the potential 



(1) 



for the atoms, and similarly with / — > b for the bosonic 
molecules. 

In the absence of interactions between particles, the 
trap levels have the energies 



(2) 



£ n = I n + - ) hui ho - 



where the principal quantum number n is positive or zero. 
In order to deal with the high degree of degeneracy of 
this potential, it is convenient to introduce the (integer) 
radial and angular quantum numbers n r and I, which are 
positive or zero, with |3l| 



n = 2n r + I. 



(3) 



Each pair (n r ,l) corresponds to a radial wave function 
R n i(r) and hence (21 + 1) common eigenfunctions of 
V f (r), L 2 and L z , 

<f>n,l,m( r ) =R nAr)Ylm(9,<t>)- (4) 

Taking into account the magnetic quantum number 

-l<m<+l, (5) 
the degeneracy of each level £ n is therefore 



fl n = i( n + l)(n + 2), 



(6) 



and the total number of states up to the shell np corre- 
sponding to the Fermi energy 



8 f 



n F ^ 

N np = 2 x n n = -(n F + l)(n F + 2)(n F + 3), 

where the factor of 2 accounts for the two hyperfine spin 
states of the atoms. In the following, it will be necessary 
to include the attractive interaction responsible for the 
pairing between fermions. Since as we show later on this 
interaction splits the degeneracy of the trap levels, we 
keep the angular momentum index in the labelling of the 
atomic energies, £ n j. The fermionic atoms are therefore 
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described by the annihilation operator c n .^ mt<T , where a 
labels the hyperfine spin state of the atom with single- 
particle wave function 4> n ,i,m(^) and eigenenergy E n ,i- 

Assuming that the atom-molecule photoassociation en- 
ergy is smaller than the trap energy spacing, ?iwho, and in 
addition that the system is sufficiently dilute that the at- 
tractive interaction between fermions is likewise less than 
Tk^ho) it is possible to tune the frequency of the photoas- 
sociation laser so as to only couple fermions in the shell 
np with Fermi energy Bp to molecules in the ground state 
of the harmonic trap. We can then ignore all shells other 
than the n^-shell for the fermions, and all trap states 
above the ground state for the molecules, which are then 
described in terms of the ground state bosonic annihila- 
tion operator b with single particle energy £q. 

Both atomic pairing and photoassociation involve the 
creation and annihilation of pairs of atoms, hence it is 
convenient to introduce pseudo-spin operators Si |20l | for 
atoms of angular momentum I in the n^-shell as 



while the number of molecules is < rib < TV. 

The pseudo-spin operator Sf 
pairs, hence any state \v) = \vi,vi 
fermions and zero molecules clearly satisfies 



annihilates atoms in 
, • ■ • , vi») of unpaired 



0. 



with 



m \v) 



n b \v) 



n\v) 



0, 



(14) 



(15) 



see Eq. <|11|) . The number operator of bosonic molecules 
has been defined by h b — b^b. V\ is the number of un- 
paired fermions of angular momentum I in each level and 
is referred to as seniority [30| in nuclear physics. 

With this formal development at hand, we now turn 
to the discussion of the fermionic pairing and of the pho- 
toassociation of atoms into molecules. It is described by 
the effective Hamiltonian 



H = (h6 + £ )tfb- 



E 



Ima 



+ V p + V am , (16) 



sr = 



Ec- 1 



-i 



i 



\l- 



Sf = - 



E "''< ■"' r ''"' • 

m— — I 
1 



(7) 
(8) 



= -(fti-nO, (9) 



where fi; =21 + 1 and hi is number operator in each level 
I. (Here and hereafter we have omitted the label np from 
the fermionic operator c nF ;i, m ,a-) They are easily seen to 
obey the SU(2) algebra 



[St i s v 



2S?S, 



i °i,v 



±sfs u , 



(10) 



Since we need only consider atoms in the Fermi level, 
the total number of relevant particles in the system is 



N 



n b 



M 



where 



M 



n b 



(11) 



(12) 



is the number of molecules (n b ) and atomic pairs (n p ) 
in the ni?-shell, or loosely speaking the number of pairs, 
and v is the number of unpaired atoms in the n^-shell. 
In that reduced Hilbert space, a complete set of states is 
given by 

|nj,n;/, • • ■ ,ni'i,n b ;v) 



1 



where M is a normalization constant. 

For a given angular momentum I, the possible number 
of atomic pairs ni in the Fermi level is 

< m < 9.1 = 21 + 1, 



where 6 is the two-photon detuning between the Ra- 
man lasers and the internal energy difference between 
atomic pairs and molecules, V p describes atomic pairing 
and V am accounts for the photoassociation of atoms into 
molecules. Before discussing these two interaction Hamil- 
tonians in detail, we first evaluate the mean-field lifting 
of the single-particle energy degeneracy, £ np — > £ nF ,i- 

In the s-wave scattering approximation, valid at T = 0, 
atoms of opposite spin interact via the two-body interac- 
tion potential 



V(r 1 -v 2 ) 



4irh 2 a 



m f 



: <5(ri - r 2 ), 



(17) 



where a < is the scattering length, negative for attrac- 
tive interactions. In the Thomas-Fermi limit, this results 
in the atoms being subjected to the mean-field potential 



V(r) = 



2nti 2 a 
m f 



where the density p(r) is given by 

p(r)^ Po (l~r 2 /R 2 TF )^, 

for 



(18) 



(19) 



r < Rtf = a\ 1Q \/2np + 3 



and is zero otherwise. Here a^o = y/^J^nf^ho is the 
oscillator length and p = {2n F + 3) 3/2 /37r 2 a^ . The 
resulting mean-field ener gy s plitting of the /-states within 
the np manifold is then 32] 



drr 2 V{r)\R nF j\ 2 

2 a 
37T a ho 



(20) 



{2n F + 3f/ 2 hu ho 



4 

3^ 



1 1(1 + 1) 
47r n 2 F 
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where we have used the WKB limit of the radial harmonic 
oscillator wave function which is valid for np 1 . For an 
atomic system to be dilute, the mean- field shift Eq. J2DI> 
should be less than the unperturbed energy £ n . This 
implies that 



1/2 



a 

— < 1, 

0-ho 



(21) 



which is equivalent to the familiar diluteness condition 
Po\a\ 3 < 1. 

It is known that for attractive short range interactions, 
the potential V(ti ~ r 2) favors the creation of a time re- 
versal state and lets pairing take place. The Hamiltonian 
V p describes this pairing correlation by including only the 
correlations for time reversal pair states 



\L = 0, M = 0:11) = (lm,l-m\00)\l,m)\l,-m), 

(22) 



m=—l 



where 



(Im, I - m\00) = {-l) l -' m (2l + l)" 1 / 2 



(23) 



is a Clebsch-Gordan coefficient. Assuming for simplicity 
that the radial part g of the pairing interaction is inde- 
pendent of the angular momenta I and I' of the atomic 
pairs involved in the interaction, V p reads explicitly 



= -*££(-i) 

l,m Z',m' 

(24) 

where the terms (— l)'~ m and (—1)' ~ m come from the 
Clebsch-Gordan coefficients (I23[l and hence account for 
the angular part of the wave functions. Strictly speaking, 
the coupling constant gi^> should be determined by the 
spatial integral 



l-m+l'-m' t J _„ ,,„,, 



4:TTtl 2 \ 



91,1' 



drr 2 Ri F Ar)Ri F ,,{r). 



(25) 



However for our present purpose it is sufficient to esti- 
mate g by replacing the spatial integral in Eq. (|25l) by 
i r n )~ 3 ^ 2 j where for a pure harmonic oscillator state in 



the np-shell the mean square radius 



is given by 



(2np + 3)a 2 from the virial theorem. This gives 



\-kTi\o\ 
m f 



(2n F + 3y 3/2 a 



-3 \ a \ -3/2 



ho 



Oh, 



tuv ho . (26) 



In this model, the strongest pairing occurs for degen- 
erate energies, £i = £y for all 1,1'. We will see that all 
fermions in the n^-shell are then coherently paired, the 
pairing energy being proportional to the degeneracy fac- 
tor il nF - In the dilute gas limit, which is equivalent to 
gil nF ~ gnp -c hoJh n, th e pairing takes then place on a 
single shell, see Ref. |32| for more details. 

Turning finally to the photoassociation of atomic pairs 
into molecules and the reverse process of photodissoci- 
ation of molecules back into atoms, we note that it is 



possible to neglect all processes involving atoms other 
than those in the n^-shell provided that their char- 
acteristic frequency, the product xV ( n b) of the pho- 
toassociation coupling constant, \i and the square root 
of the mean number of molecules, remains small com- 
pared to the frequency separation between neighboring 
shells of the trap. This condition is well fulfilled for 
typical laser strengths and trap depths, in which case 
xV ( n b) — 10 2 n f 1 ^ 2 s _1 -c ^ho — 10 3 s _1 . In this case, 
the atom-molecule coupling can be approximated by the 
Hamiltonian 



V a 



= £<-*: 



l,m 



where we only include the coupling between time-reversal 
atomic pair and bosonic molecules in the ground state 
of the harmonic trap. The photoassociation coupling 
constant \ is proportional to the far off-resonant two- 
photon Rabi frequency associated with two nearly co- 
propagating lasers with frequencies u)\ and u>2 Q, 



X 



Xo 



47T 

Xo 



drr 2 R { Q b l(r)Rl F ,(r) 



3/2 3/2 ; 



(28) 



where i? b Q is a ground-state wave function of molecules 
of mass of 2m/ with a spatial width about aho, and in 
the second equality we have again replaced the radial 
wave function of the atoms by (?" 2 F ) _3 / 2 - From the cou- 
pling constant between 87 Rb2 molecules and 87 Rb atoms 
of Ref. (3j| we estimate the photoassociation coupling 
constant to be of the order of xo = 7.6 x 1CP 7 m 3 / 2 s _1 . 
In terms of the pseudo-spin operators, and with Eqs. 
and (J23, the total model Hamiltonian Eq. lO fi- 
nally reads 

h = {n8 + £ )bh+Y,^n F ,i{S! + n l /2) 
i 

+ x £(S+6 + 6tsp_ fl £s+s-. (29) 
i i,i' 

This Hamiltonian clearly conserves the spin operators Sf. 
Applying this operator on the state \v) we have 

S» = Si{S t + l)\v) 

= {s+Sf+sf(sf-i)}W) 

= {Qi/2- vi/2)(ni/2-vi/2 + l)\u), 

where we have used the identity S 2 = S^Sf + Sf (Sf — 1) 
and Eqs. (JUJl and fT5|t. so that 

Si = -(ili - vi). 
This allows us to identify the operator 



Si 



(30) 
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as the operator giving the total number of fermionic pairs 
and molecules, or loosely speaking the "pair number" op- 
erator, which is easily seen to also be a conserved quan- 
tity. 

In the limit g — > the Hamiltonian (|29|) reduces to 
the Tavis-Cummings Model pll 12*2^ of quantum optics, 
while for \ ~ * 0, it becomes the Pairing Model 0, 
l35| . for which Richardson first gave an exact solution 
in the context of nuclear physics. We also mention a 
recent family of exactly solvable models of atom-molecule 
proposed in Ref. [3(1 . 

The discussion of the following sections concentrates 
specifically in that situation where the mean-field en- 
ergy shift in single-particle energies is smaller than the 
photoassociation coupling. In this degenerate model, the 
single-particle energies of all atoms in the n^-shell can 
then be taken to be equal, £ nF ,i = £f, and the Hamilto- 
nian (|29(l simplifies to 

H = -huj (S z + S) + U X (S+b + tfS-)-hgS + S-, (31) 

where we have introduced the total spin operator 



S = E S *> 



(32) 



oj = 5 + (£ - 2£ F )/K, {l np = J2i &h v — v h so tnat 
the total spin is 



s = J2 s < = Un nF -v) 



(33) 



In Eq. (|3 1|) . we have neglected constant terms propor- 
tional to the conserved quantity M. 

Assuming an oscillator length aho = 3.2 x 1CP 6 m for 
the mass of 6 Li, a = —114 nm for its scattering length, 
and Who = 1000 s~\ then xn 3 J 2 ^ 10 2 s -1 . The validity 
of the degenerate model, \£ nF ,n F — £n F ,l=o\ < X n F re- 
quires np < 10, which corresponds to a total number of 
fermions N n „ < 10 3 . 



A. Pairing Model 

For a total number of particles N = 2M + v < O , 
which corresponds to the Fermi energy shell less than 
half-filled since there are two hyperfine atomic states in- 
volved, the seniority v can take the values 



0,2,4, ...,N {N even) 
1,3,5,..., N (N odd), 



(35) 



while for O < N < 20, corresponding to a shell more 
than half filled, the permissible values of v are 



0,2,4, ...,20- N (N even) 
1,3, 5,..., 20 -TV (N odd). 



(36) 



The degenerate Pairing Model of Refs. [33, |3l| is 
obtained by neglecting the photoassociation coupling, 
X = 0, in the Hamiltonian (|31() . In this case, the to- 
tal energy is given as a function of seniority v and the 
total number N of fermions by 



E, 



f(AT 



v)(2Q-y-N + 2), (37) 



and it is minimized for v — or 1. 

For an even particle number, the ground state corre- 
sponds therefore to all pseudo-spins aligned, S — 0/2, 



E pm (N, v = Q) = -loN - gN(2tt - N + 2). 



(38) 



The first excited state corresponds to v — 2 unpaired 
atoms, or S = 0/2 — 1, and its energy is 



E pm (N, u = 2)- E pm (N, v = 0) = gSl. 



(39) 



Thus, the energy needed to break up an atomic pair into 
two un-paired fermions is independent of the number of 
fermions in the n^-shell. This energy is consistent to Bo- 
goliubov quasi-particle energy based on BCS variational 
state except for corrections of relative order I/O. 



III. GROUND STATE 

In this section, we discuss the dependence of the 
ground state of the model Hamiltonian H31(l on the ratio 
of number of paired fermions and molecules, M = n p +rib, 
to the degeneracy O of the Fermi level (we drop the sub- 
script 'hf' for notational clarity from now on). We iden- 
tify qualitatively different types of ground states, a pair- 
dominated ground state and a molecular-dominated one, 
MS i\ function of the parameter 



(34) 



where uj = 5 + (£q — 2£p)/h is the photoassociation fre- 
quency detuning. 



B. Photoassociation 

In the presence of photoassociation, x / 0, the eigen- 
states of the system consist of a coherent superposition of 
atoms and molecules. For fixed N, we can classify them 
by the total spin S = (O — v)/2 : each manifold consisting 
of (M + 1) eigenstates, where M — (N — v)/2 — n p + rib 
is the total number of pairs. 

For positive detuning, u> > 0, the ground state is a pure 
fermionic state, i.e., N = 2n p + v. On the other hand, 
for a large negative detuning energy, uj < 0, and even 
particle number, the ground state is reached when all 
particles are molecules, corresponding to a zero seniority, 
or maximum spin state. 

Figure^shows the energies E v of a few seniority states 
relative to that of the state v = as a function of 
k = oj/xVO- This example is for r\ = gy/Q/x — 1 and 
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FIG. 1: Excitation energies, in units of XV^, relative to 
the ground-state energy for the seniority states v = 2 (solid 
line), v — 4 (dashed line), and v — 6 (dash-dotted line), as 
a function of the dimensionless parameter re = w/xV" for 
■q = gVtl/x = 1 and N = Q = 120. 



a half-filled Fermi energy shell, N = ft = 120, corre- 
sponding to rip = 14, and results from a direct numerical 
diagonalization of the Hamiltonian. It shows that over 
the wide range of k considered here, the ground state 
is always the maximum spin manifold v = 0. For large 
negative detunings, the energy differences E v — Eq are 
v jl times the single-particle energy difference between a 
molecule and two fermions, and increase with \uj\. For 
large positive detunings, on the other hand, the energies 
of the excited states approach the values v/2 times gft, 
see Eq. (I39f) . In the crossover region where the nature of 
the ground state changes from atomic to molecular, the 
ground state is a coherent superposition of atoms and 
molecules that is likewise the maximum spin state. We 
have verified that in the absence of pairing interaction, 
g = 0, the ground state is likewise the state of maximum 
spin (not shown in figure). 

In the following we concentrate on the state of max- 
imum total spin, \S = ft/2), for an even number of 
fermions, N = 2M. The eigenstates of the atom- 
molecule system in the spin manifold S = ft/ 2 have the 
general form 



A I 



= ]T C x (n p )\S = fl/2,S z = -S + n p ) f 

n b =0 

(g)|n h = M - n p ) b , 



(40) 



where n p denotes 



a number of fermionic pairs and 
A = 0, 1, • • • M represent eigenmodes of the system with 
eigenenergies E\. We have found these eigenstates and 
the associated eigenenergies by direct diagonalization of 
the Hamiltonian (|31l) . 

The average number of molecules in the ground state 
is shown in Fig. [3 as a function of k = uj/xVfi and M 
for g = and ft = 120. As expected, for large positive 



FIG. 2: Normalized average number of molecules in the 
ground state as a function of re and M for fl — 120 and 
5 = 0. 




FIG. 3: Normalized joint coherence function, \G ar , 
functions of re and M for fl — 120 and for g — 0. 



detuning ui > 0, the ground state population consists al- 
most exclusively of atoms, while it is mostly made up of 
molecules for large negative detunings. In the crossover 
region around u = 0, the ground state consists of a co- 
herent superposition state between molecules and atomic 
pairs. 

It is possible to characterize this superposition in terms 
of the normalized joint coherence function 



Gar. 



2(S+b) 



MVft 



(41) 



which is shown in Fig. as a function of k and the num- 
ber of pairs M. The joint coherence of the atomic and 
molecular fields shows a remarkable enhancement in the 
crossover region as well as a change in shape as a func- 
tion of M, due to the nonlinearity of the atom-molecule 
coupling (|27fl . This dependence on the filling factor M/ft 
can be understood more quantitatively by considering the 
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two limiting cases M/Q <C 1 and M/f2 ~ 1. 



2. M/Q ~ 1 — Mapping on a binary atomic-molecular 
condensate 



1. M/Q <C 1 — mapping on a linear coupled-boson system 

In the limit of small filling factors, it is convenient to 
describe the system in terms of the Holstein-Primakoff 
mapping [3^| of the SU(2) generators S + = (S~Y and 
S z in terms of bosonic operators. According to this map- 
ping, the Hilbert space of the group SU (2) is carried by a 
subspace of the bosonic Fock space given by the bosonic 
vacuum |0)d and the bosonic operators d and <v, with 

[d,rf] = l, d|0> d = 0, (42) 

The bosonic space being spanned by the n^— boson states 

1 



\rid)d 



■A<$T d \G)d fern* = 0,1, 2,. 



,M. 



Since we restrict our considerations to the subspace char- 
acterized by the angular momentum quantum number 
S = 0/2, we can map the operators S and S z as 



s+ -> VfidU i 



dU 

~n~' 



s~ 
s z 



rffd 

7T 



(43) 



In the limit M/Q < 1, we only need a lowest-order of 
the operators @3J), that is, replace the square roots by 1. 
The Hamiltonian reduces then to that of a linearly 
coupled two-mode boson system, 



H —> H] 



linear 



= -{L0 + gQ)d)d + X ^(d)b + tfd), (44) 



with the constraint M = rid + rib. 

This Hamiltonian can be diagonalized by the transfor- 
mation 

c_ = cos 9d^ — sin 9b^ 
= sin Q<$ + cos 9b\ 

with cot 29 = (w + gQ) /2xVQ, ,0 < 9 < 2tt, to give 



linear 



(45) 



with energies 



- (cj + gfi) cos 2 - 2x VO cos sin 0, 

- (oj + gQ) sin 2 6* + 2xVQ cos sin 0. 



Since e_ < e + for \ > 0, the ground state is 
(l/A/M!)(cL) M |0) d |0) b . 



Tikhonenkov and Vardi j26j] showed for in the homoge- 
neous case that when the total number of pairs is equal 
to the available momentum states of the fermions, the 
system of fermionic atoms and bosonic molecules can 
be mapp ed onto a two-mode atomic-molecular BEC sys- 
tem |4Ct lilt l42j . The corresponding situation in our case 
occurs for M — Q, with an additional atom-molecule 
two-body collision term required in addition. 

To show how this works we introduce the two-mode 
Hamiltonian of a two-component condensate of atoms 
and molecules, 

H am = -ujblb m + ^(bib a ba + h.c.)-hlb m bib a , (46) 

where b\ and 0^ are the bosonic creation operators for 
the atomic and molecular modes, respectively. Clearly 
the total number of particles N am = n a + 2n rn is con- 
served, where n a and the number of atoms and 
molecules, respectively. For N am even, a general state of 
the system can be expressed as 

M 



\$am) = ^2 C{m a )\n a — 2m a ) ® \n m = M — m a ) 
ro o =0 
M 

= C{m a )\2m a ,n m ), (47) 

m a =Q 

with M = N am /2. In this representation, the matrix 
form of the Hamiltonian l|46(l is 

(2m a ;n m \H am \2m 

= -u>(M - m a ) - gm a (M - m a ) 
(2(m a - l);n m + l\H am \2m a ; 

= (2m a ; n rn \H am \2(m a - l);n m + 1) 

= xV m a(m a - l/2)(M - m a ). 

Similarly, for our model the matrix form of the Hamilto- 
nian iPfljl in the "pair number" representation is 

(n b ;n p \H\n b ;n p ) 

= -u(M - n b ) ~ g(n b + AM + 1)(M - n p ) 

(n b - l;n p + l\H\n p ; n b ) 

= (n p ; n b \H\n b -l;n p + 1) 

= xVM^b + AM + 1)(M - n b ). 



where 



AM = Q- M, 



(48) 



with AM > 0. In the limit M ~ Q, that is, when we can 
neglect AM/M <§; 1 and other terms of order of 1/Q, 
these two Hamiltonians are same under the transforma- 
tions m a «-> n b and n m <-> n p . 
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FIG. 4: Lowest excitation energies as a function of k for 
M/Q, = 1 (solid line) and M/Q = 0.5 (dashed line), and 
M/Q = 0.1 (dash-dotted line), at fixed pair number M = 120, 
where g = 0. The asterisks correspond to a linear approxima- 
tion for small M/Q (see text). 



3. Intermediate regime 

The behavior of the system apart from these two lim- 
iting cases deviates from both models. The most striking 
difference between these regimes appears in the energy, 
A, of the first excited state. Figure 0] shows A/xVQ in 
the absence of pairing interaction for M/Q = 1 (solid 
line), M/Q = 0.5 (dashed line), and M/Q = 0.1 (dash- 
dotted line) as a function of k = lu/xVQ for fixed 
M = 120. In the case of M/Q = 0.1, A agrees well with 
the one-particle energy difference e + — e_ = lo 2 + A\ 2 Q 
of the linear coupled-boson model (|45|l (shown as aster- 
isks in the figure). Increasing M/Q shifts the location 
of the minimum of the energy gap and reduces the value 
of its minimum. For M — Q, finally, the minimum gap 
approaches zero at k ~ 2, consistently with a transition 
point in the atom- molecule condensate system |4Ct 1 121 . 
This result is indicative of the appearance of a quantum 
phase transition in the limit M — > oo. 



C. The role of the pairing interaction 

We now examine in more detail the ground-state statis- 
tics of the molecular field in the presence of pairing inter- 
action V p . Since the cases of M <C Q and M ~ Q can be 
mapped onto relatively well-known systems, we present 
results for the situation of a half-filled shell, M/Q = 0.5, 
only. Figures El and show the probability P(n&) of hav- 
ing rif, molecules in the trap, the molecule statistics, for 
2M = Q = 120 as a function of the dimensionless pa- 
rameter k, with g — and rj = gy/Q/x = 10, respec- 
tively. There are several quantitative differences between 
the two cases: the transition from an atomic to a molec- 




FIG. 5: Molecule statistics P(nt) as a function of the 
dimensionless parameter k in the absence of pairing in- 
teraction, g = 0, for 2M = Q = 120. 




FIG. 6: Molecule statistics P(rib) as a function of the 
dimensionless parameter k in the presence of pairing in- 
teraction, r) = gV&~/x = 10, for 2M = Q= 120. 

ular ground state is shifted by the pairing interaction, 
and the width of the crossover region in detuning space 
is significantly broader. 

This behavior can be understood by noting that the 
pairing interaction gives rise to an additional detuning 
effect depending on the number of molecules, as shown 
by the diagonal terms of the Hamiltonian (|48|l with 
M = AM = Q/2. This interaction leads to number- 
dependent energy shifts and dephasing between different 
molecular number states. A similar effect has been stud- 
ied in the context of a Jaynes-Cummings-like description 
of photoassociation and has been referred as "nonlinear 
detuning" Q . 

We can estimate the position of the "resonance" point 
^res(g) where the ground state goes from being molecular 
to atomic in nature by taking (rib = M/2+1; n p — M /2 — 
l\H\n b = M/2 + l;n p = M/2 - 1) = (n b = M/2;n p = 
M/2\H\n b = M/2;n p = M/2). This gives 

Q 

Wres(g) = -g-^, 
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where we have neglected a term of order 1/17. 

The shift in uj ves (g) due to the pairing interaction is 
further illustrated in Fig. which shows the magnitude 
of the joint coherence function, |G om |, as a function of 
k = lu/xVQ for three values of the pairing coefficient g. 
Finally, Fig. [S] shows the entanglement entropy E{pb) 
of the ground state, obtained from the von Neumann 
entropy of the molecular reduced density operator |43| 

Pb = Tr f (p), 

as 

M 

E(Pb) = -J2 \C(n b )\ 2 log\C(n b )\ 2 , (49) 

n b —0 

where the logarithm is taken in base 2, for three values 
of the pairing interaction strengths, rj = 0, 5, 10. The 
entanglement in Figure |5J is divided by a maximum en- 
tanglement, logM. Consistently with the results of of 
Figs. [5] and for the molecular statistics, the pairing 
interaction reduces the entanglement. 

IV. DYNAMICS 

In the absence of pairing interaction our model is 
equivalent to the Tavis-Cummings model, whose dynam- 
ics has been studied in detail in the context of coherent 
spontaneous emission from a system of N two-level atoms 
interacting with a quantized radiation field |U Q ■ The 
main purpose of this section is to extend this work to the 
study of the system dynamics in the presence of V p . One 
important result is that the nonlinearity of this interac- 
tion produces a self-trapping transition. 

We assume that the system consists initially either of 
atomic pairs or of molecules only, corresponding to a 
maximal spin state. A photoassociation beam is applied 
from t = on, and we examine the subsequent coher- 
ent dynamics of the system. Note that although this 
problem resembles the dynamics of a bosonic Josephson 
Junction [2ll |2S| in an asymmetric trap for the initial 
imbalance in populations, the nonlinear coupling term in 
our model leads to considerably richer dynamics. Since 
the total spin is a constant of motion of the Hamilto- 
nian i|31|) , we confine our discussion to the maximal spin 
state \S = 17/2). 

A. Semiclassical approximation 

Before proceeding with a full quantum analysis, we 
first consider the semiclassical approximation. Our ap- 
proach is very similar to that taken in Ref. |4l|. Intro- 
ducing the two operators 

J+ = S + b, 

J- = tfS-, (50) 



0.9 
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FIG. 7: Normalized joint coherence, |C? a m|, of the ground 
state as a function of k for three values of the pairing 
interaction strengths rj — gVTi/x, indicated in the insert. 
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FIG. 8: Entanglement entropy, E(pt), of the ground 
state, in units of the maximum entanglement logM, as 
a function of k for three values of the pairing interaction 
strengths r\ — g\/Q./x indicated in the insert. 

results in the Heisenberg equations of motion 

i^-J+ = ujJ + + 2 X b'<bS z + X S + S- -2g,J+S z , (51) 
at 

ij t J- = -ujJ- -2 x b f bS z - xS + S- + 2 5 S' Z J_,(52) 

ij f S z = X(J+ - J-)- (53) 

Here we note the following relations between operators 
and conserved quantities, 

S+S- = S{S + 1) - S Z (S Z - 1), 6 f 6 = M — S — S z . 

As usual we introduce a semiclassical approximation 
by factorizing the mean values of the various opera- 
tors that appear in the Heisenberg equations of motion, 
such as (S Z J+) = (S Z )(J+) and (J-S z ) = (J-)(S Z ), 
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etc. Introducing the c— number functions s z = (S z ), 
j x = ((J + +J_))/2, and j y = (( J + - J_))/2i, and neglect- 
ing corrections of order l/f2, that is, setting (S + S~) = 
S 2 — s 2 , we obtain the semiclassical equations of motion 



dt Jx 
d . 

Jt Jy 
d 

Jt Sz 



where h(s z ) 
additional conserved quantity 

d_ 

If 



= -ujx + 2gj x s z - xH s *) 
3s 2 + 2(M - S)s 



S 2 



(54) 
(55) 
(56) 
Noting the 



2x 2 X 



0. 



we find that the coupled equations l|54|l are equivalent to 
the classical Newtonian equation of motion 



dt 2 ' 



ds z 



U(s z 



(57) 



where the potential U(s z ) is determined by the initial 
conditions for j x (fi) and sq = s z (0). It is sufficient for our 
purpose to assume an initial Fock state, so that j x {0) = 
and U(s z ) is given by 



U(s z ) 



9 2 4 / 



gus Q 



g 2 (s ) 2 + 2 X 2 (M-S) 



[2x 2 S 2 - uj 2 s + togs 2 ] s z 



(58) 



Since the potential has a quartic form of s z , Eq. I|57|) 
can be solved analytically in terms of Jacobian elliptic 
functions. The derivation of the general solutions, which 
is straightforward but lengthy, is given in Appendix A. 



B. Coherent dynamics, g = 

We first examine a typical behavior of the semiclassical 
dynamics in the absence of pairing interaction, 5 = 0, for 
M/n < 1, M/n - 0.5, and M/fi ~ 1. 

From the semiclassical solutions (|A8() in the case of 
9 = and on the exact resonance uj = 0, the population 
imbalance between fcrmionic pairs and molecules, (n p ) — 
(lib) = 2S — M + 2s z , and the coherence function j y is 
given by 



(n p ) - (rib) 



M 



= 2sn 2 (xVnt + fak) - 1, (59) 
= 2 sn(xVO,t + <f>; k) X cn^VOt + <f>; k) 
x dn(xVcit + (/); k), (60) 



where sn, cn, and dn are Jacobian elliptic functions |45j |. 
Here, k — ^/M/Q is the elliptic modulus and is a 



phase factor determined by the initial conditions. It is 
equal to <j) — K, corresponding to the complete elliptic 
integral of the first kind, for an initial fermionic Fock 
state | rip = M), and to = for an initial the entire 
population being in the molecular mode. We note that 
j x (t) is zero for all time on the exact resonance. 



1. M/Q <C 1 — Linear coupled-bosons regime 

In this case, the elliptic modulus k <§; 1 (see Ap- 
pendix A) so that the elliptic functions in Eqs. I|59l) and 
(|60|l can be approximated by sn(it, k) — > sinu, cn(u, k) —> 
cosu, and dn(u, k) — > 1, respectively. The imbalance in 
atomic and molecular populations which undergoes Rabi 
oscillations at the frequency 2xVQ is given by 



({rip) - [n b ))/M = cos (2 X Vnt) 
for an initial fcrmionic pair state, 



((n p ) - (n b ))/M = -cos(2xVm) 

for an initial molecular state. These results are equivalent 
to those obtained directly from the Hamiltonian (|44|l . Wc 
have compared these solutions with the exact quantum 
mechanical dynamics obtained numerically, and checked 
that the linear approximation agrees with the quantum 
results for M/Q, < 0.2 and for times shorter than t ~ 

Tr/xVn. 



2. M/Sl ~ 0.5 — Intermediate regime 

Figure shows the normalized population difference 
((n p ) — (h},))/M, and Fig. |5Jd shows the normalized co- 
herence function 2j y /M \j£l/2, as a function of the di- 
mensionless time r = xV&t for a system initially either 
in a pure atomic state or a pure molecular state and for 
M/f2 = 0.5. The circles correspond to the semiclassi- 
cal description, while the lines are the results of a full 
quantum-mechanical analysis. The anharmonicity due 
to the nonlinear atom-molecule coupling is clearly ap- 
parent, and also shows that the semiclassical dynamics 
approximate the quantum dynamics very well. 

We note that the atomic pair state in the half-filled 
shell corresponds to a Dicke superradiant state [2^ |. 
which is known from quantum optics to give rise to the 
strongest collective enhancement of transition probabil- 
ities. This enhancement is proportional to the prod- 
uct of the number of particle pairs M and the number 
of hole pairs Q — M, and is maximum for M = Q/2. 
From Eq. Q59[l . we have that for sufficient short times 
t = t/xVfi < 1 the average number of molecules builds 
up as 

(n b ) = cn 2 (xVQt + 4>) 

^ T «i M (Q - M) { x tf = ^n 2 t 2 . 
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FIG. 9: Comparison of the semiclassical and quantum dy- 
namics for 2M = O = 120, U) — 0, and g — 0. Figure (a) 
shows the population difference between fermionic pairs and 
molecules {{n p ) — {hb))/M. Figure (b) plots the coherence 
function 2j y / M\/~Q. as a function of the dimensionless time 
r = xV&t. The solid and dashed lines give the quantum re- 
sults for an initial fermionic pair sate and a pure molecular 
state, respectively. The corresponding semiclassical results 
are indicated by filled and open circles, respectively. 



3. M/Q. ~ 1 — Binary atomic-molecular BEC 

As shown in Ref. [2(J, the coherent dynamics in this 
regime is qualitatively very similar to that of binary con- 
densate of atoms and molecules. For M = Q, corre- 
sponding to k — 1, the population imbalance between 
atomic pairs and molecules is given in the semiclassical 
approximation by 



- (n b ) 
M 



= 2 tanh 2 (xV^t 



-1, 



(61) 



indicating that the point (n p ) = M is stationary. How- 
ever, the system is dynamically unstable against small 
fluctuations , see Ref. f° r a detailed discussion in the 
context of binary condensates of atoms and molecules 
and Ref. ^(| in the context of second-harmonic genera- 
tion. 



C. Self-trapping transition and quantum dynamics 

To conclude, we discuss the coherent dynamics of the 
system in the presence of pairing interaction, g ^ 0, con- 
sidering only the case of half-filling for simplicity. We 
consider specifically the example 2M = ft = 120, which 
corresponds to nj? = 14, and take an initial state as the 
Fock state \n p = M;rif, = 0). Figure ITTI of the appendix 
shows the phase diagram of the semiclassical dynamics 
in the k — nsp&ce, illustrating that a self-trapping tran- 
sition |2?l [2|j takes place when varying the dimension- 
less detuning frequency k = uj/xV^: provided that the 
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FIG. 10: Schematic potential curves, U{s z ) — U{sq), around 
a self-trapping transition point kq. These curves are at (a) 
above ko, (b) at transition point k = /to, and (c) below kq. 
The filled circle indicates an initial position of classical parti- 
cle. 



pairing interaction strength r\ = g\ZTi/x exceeds a criti- 
cal strength r\ c . For the case of half-filling case, we find 
ry c = 5.0302. 

This transition can be interpreted physically from the 
motion of "classical particle" in the "potential" l|58|) . Fig- 
ure EH displays schematic potential curves, U(s z ) — U(sq) 
as a function of s z in the vicinity of a self-trapping transi- 
tion point kq- For our specific initial conditions the par- 
ticle "velocity" is initially zero, ds z (0)/dt oc j y (0) = 0. 
For k > kq, Fig. 1101-fa). the classical particle oscillates 
periodically in the potential. As n approaches kq, one 
additional potential barrier appears, and at the transi- 
tion point its height equals the initial potential energy of 
the particle (see Fig. HUI-Cb)). At that point, the particle 
rests on the potential maximum after reaching it. Below 
the critical point, the barrier confines the particle in a 
narrow range, as shown in Fig. 1101 -fc). The self-trapping 
effect provides a sudden suppression of the amplitude of 
coherent oscillations. Since the key factor in achieving 
this transition is the quartic term in the potential JSSJ, 
it disappears in the absence of pairing interaction. 

Fig. II II shows the time evolution of the population dif- 
ference {{hp) — {hb))/M for r\ — 6 and for several detun- 
ing energies. The semiclassical solutions (dashed lines) 
clearly show the self-trapping as a sudden suppression of 
coherent oscillations from just above to just below the 
transition detuning. The dotted line in Fig. HH-fc) cor- 
responds to the semiclassical solution for the threshold 
detuning kq = —3.2307. The solid lines show the exact 
quantum solutions. Apart from the transition point, the 
quantum and semiclassical dynamics are similar at least 
for short enough times. However, the oscillations of the 
quantum solution deviate from those of the semiclassical 
solution near the transition point, Fig. Illf b-c). Since in 
the semiclassical picture the height of the potential bar- 
rier near the transition point is just below or above the 
initial potential energy of a particle, the quantum mo- 
tion of the particle is very sensitive to fluctuations and 
hence deviates significantly from its classical counterpart. 
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FIG. 11: Population imbalance versus the rescaled time r = XV^t f° r V = 6. Quantum (solid line) and semiclassical (dashed 
line) solutions are shown, respectively, for the detuning parameters (a) K = and (b) k — —3.1, and (c) re = —3.3, and (d) 
k = —4.0. The dotted line in (c) corresponds to the semiclassical solution at the transition point Ko = —3.2307. 



The initial Fock state provides fluctuations of coherence, 
and large quantum fluctuations of the population imbal- 
ance arise as a result. We have verified numerically that 
the number fluctuations near the transition point are en- 
hanced by an order of magnitudes as compared to those 
far away from that point. 



Figs.ll2landll3lcompare the quantum and semiclassical 
time-averaged population imbalance as a function of re for 
r\ = 2.0, and for r\ = 6.0, respectively. In contrast to the 
second case, there is no self-trapping transition in the first 
case. Hence the semiclassical time-averaged population 
imbalance (cross) is a smooth function of re, and agrees 
well with the quantum results. For the strong pairing 
coupling of Fig. ^| in contrast, an abrupt jump of the 
semiclassical time-averaged value occurs when varying re, 
a signature of the self-trapping transition. Due to the 
large quantum fluctuations, it differs markedly from the 
time-averaged quantum result near a transition point. 



V. SUMMARY 



In this paper, we have considered the coherent pho- 
toassociation of fermionic atoms into bosonic molecules 
trapped in a spherically symmetric harmonic trap. We 
showed that under a realistic set of conditions this sys- 
tem can be mapped onto a Tavis-Cummings Hamiltonian 
with an additional paring interaction using pseudo-spin 
operators. We carried out an exact numerical diagonal- 
ization of the Hamiltonian to determine the ground state 
of the system, investigating the crossover from a predom- 
inantly atomic to a predominantly molecular state. We 
also investigated the joint coherence and the quantum en- 
tanglement between the atomic and molecular fields, and 
found that the atomic pairing interaction suppresses the 
entanglement between fermions and bosons. We then an- 
alyzed the coherent dynamics of photoassociation due to 
the nonlinear atom-molecule coupling. Using a semiclas- 
sical factorization ansatz, we showed the appearance of a 
self-trapping transition in the presence of pairing interac- 
tion. An exact quantum solution illustrated the impor- 
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FIG. 12: Time- averaged population imbalance as a func- 
tion of K for r\ = 2.0. Quantum result (cross) and semi- 
classical result (open circle). 



FIG. 13: Time-averaged population imbalance as a func- 
tion of k for t) = 6.0. Quantum result (cross) and semi- 
classical result (open circle). 



tant role of quantum fluctuations in the neighborhood of 
that transition point. Future work will extend this study 
to a detailed analysis of the non-degenerate model and to 
multi-well superradiant systems. For instance, preparing 
an atomic Fermi gas in a Josephson-type configuration 
and applying a photoassociation beam should lead to the 
efficient production of spatially correlated molecules. 



state as an atomic state in the maximum spin manifold 
S = 0/2. The extension to other initial states is straight- 
forward. 

1. Case g = 

In this case, the "potential" U(s z ) is a cubic function 

of s z , 
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APPENDIX A: ANALYTIC SOLUTIONS IN 
TERMS OF JACOBIAN ELLIPTIC FUNCTIONS 

In this appendix, we obtain analytic solutions of semi- 
classical dynamics that obey the Newtonian equation of 
motion Eq. (|57|1 . and show the phase diagram of the semi- 
classical dynamics based on those solutions. 

The solution of Eq. I|57[) has the general form 



t = 



s y/2[U{s ) - U{s z )\ 



(Al) 



with the potential U{s z ) given by Eq. 1)58(1 . 

We analyze the solution for the two cases g = and 
g =/= 0, separately, and for simplicity we take the initial 



J 



U{s z ) = -2 X 's% + 



■+2 X 2 s a s 2 z +[2 X 2 S 2 -uj 2 s ]s z . 

(A2) 

Here we have used the initial condition sq = —S + M. 
By introducing the normalized quantities 

s z = s z /2S (-1/2 <s z < 1/2), s = s /2S, 

we obtain the explicit form of the denominator of the 
right-hand side of Eq. (|A1|) , 

<p{a z ) = 2[U{s )-U{s z )} 

= A x 2 {2Sf{s z - s )(s z - s+){s z - s_) (A3) 



K 1 



s± = — ±-Vk 4 -16s q k 2 + 16, 

o o 



(A4) 



where k — uj/xv2S. We note that the variables s± 
are always real- valued for any magnitudes of k 2 , because 
-1/2 < s < 1/2, and then s + > s > s z {t) > s_. The 
solution can be obtained by integrating the form 



ds 



vVO*) xV 2S ( s + 




(A5) 
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where, 



4>{t) = arcsin \ 



(s+ - s_)(s - s z ) 
(s - s_)(s+ - s z ) ' 



The integral that appears in that equation is an ellip- 
tic integral of the first kind. Noting that the integration 
within < < tt/2 gives rise to a complete elliptic inte- 
gral of the first kind, K, we find 



X\/25(s + -s-)t + K = sn^ 1 (sin 0(f) ; fc) (A6) 

where the function sn^ 1 is the inverse of the Jacobian 
elliptic function and k — y/(s~o — s -)/(&+ — s_) denotes 
the elliptic modulus. This gives the evolution of s z (t) 

s z (t) = s- + (s - s_)sn 2 (xV 2 S(s+ ~ s-)t + K; fc) . 

(A7) 

At the exact resonance, lu = 0, we have that s+ = 1/2, 
S- = —1/2, and k = \J M/2S so that this expression 
reduces to 



s z (t) = -S + Msn 2 ( X V2St + K;k). (A8) 

In terms of s z (t) the coherence functions j x (t) and j y (t) 
are given by 

. . lu , . , . 1 ds z (t) . . „. 

J,(*) = ^-.(t), Jv(*) = ^-^ i - (A9) 

This results in the expressions for the difference in atomic 
and molecular populations (n p ) — (rib), and the coherence 
function j y of Eqs. and (|rJU)l . respectively. 



2. Case 3 / 

The presence of pairing interaction renders the poten- 
tial quartic in s z , see Eq. (|58|) where we have again as- 



sumed that the initial state is an atomic state of maxi- 
mum spin, \S, ~S + M). It is convenient to introduce the 
function f(s z ) as 



L?(s z ) = -g 2 (2S) 4 (s z -s a )- f(s z ), 



(A10) 



where 



f(s z ) 



as 



(8* 



+ (3s z + 7 
si)(s z - s 2 )(s 2 



s 3 ), (All) 



s — 2k/t] - 4/?7 2 , f3 = —Sq + k 2 /t] 2 , and 7 



2kSq/ti — k 2 so/t] 2 + I/77 with rj = gyQ/X- 

The variables s\, §2, and S3, which correspond to the 
roots of the cubic equation, f(sj) = 0, are obtained by 
"Cardano's formula". With £ = e l27r / 3 , p = -a 2 /3 + 0, 
and q = 2a 3 /27 - a/3 /3 + 7, and also D = -(4p 3 + 27<? 2 ), 
those roots are given by 



D 
3~ 



1/3 



1/3 



(A12) 



2 6 V 3 
where j = 1, 2, 3. 

The numbers of real and complex roots are determined 
by the sign of the polynomial discriminant D. If (a) 
D > 0, all three roots are real and unequal. If (b) D < 0, 
one root is real and two are complex conjugates. If (c) 
D = 0, two roots are equal for q ^ 0, and all roots are 
equal for q = p = 0. 



a. Case D > 

Suppose that sq > s z > s a , and s a > Sb > s c or 
Sb > s c > So, where each s Qi b !C corresponds to one of the 
roots Sj's. The solution of Ea. l|Al|) reads then 



*=/ ~4^ = j , = ■ (A13) 

s V^ 5 ) gSy/{so - s b ){s a - s c ) Jo Vl-& 2 sin 2 6> 



where 



so that 



. /(s - s c )(s - s z ) (so - s a )(s b - s c ) 
<j)(t) = arcsin W — ^— — , * = W 7= rr, (A14) 

V ( s - S a )(s 2 - S c ) U (S - Sb)(Sa - Sc) 



so(s a - s c ) + s c (s - s a )sn 2 \ -gSy/(s - s b ){s a - s c )t; k) 

s z (t) = j-± —^L. (A15) 

(s a - s c ) + (s - s a )sn 2 {-gSy/(s - s b )(s a - s c )t;k\ 



Case D < s a ) {(s 2 — sr) 2 + s 2 }. With the change of variable from 



Letting s a label the real root and with Sb — s* = 
sr + is j, we have that (f(s z ) — —g 2 (2S) i (s z — Sq)(s z — 
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s z (so > s z > s a ) to as 
Sq — s z A 1 — cos (f> 



Bl + cos0' k = x / g> »> , (A17) 



( A16 ) . _ /(5o-s a ) 2 -(A-B 



where A = ^/ (s — s^) 2 + s 2 > and B = 

J (s a — Sfl) 2 + s 2 > 0, and taking the elliptic modulus the integral Eq. (|A1I) can then be replaced by 



ds 



dd 



I s vW) 2gSVAB Jo y/l - k 2 sin 2 9 
The semiclassical solution is then given by 

s a A + s B - (s a A - s B)cn (-2gSy/ABt; k 



(A18) 



A + B -{A- B)cn -2gSVABt; k 



(A19) 



c. Case D = 

In this case, the solution of Eq. IjAlfl can be expressed 
in terms of elementary functions. For q ^ the solutions 
are equivalent to Eq. I|A15(I in which the elliptic functions 
are replaced by trigonometric functions for k = (st = 
s c ), or hyperbolic functions for k — 1 (so > s a = Sfc). If 
q = p = 0, the function /(s z ) has triple degenerate roots 
at a point s a = sj, = s c , and the corresponding solution 
is obtained by 



corresponding to Eq. (IA15I) differ from Eq. 1A19H . This 



,(*) = Sa + 



1 - {-gS(sa - s a )ty 



(A20) 




3. Phase diagram of semiclassical dynamics 

In this subsection, we discuss the structure of the semi- 
classical dynamics in n — rj parameter space for the spe- 
cific case of a half-filled shell, 2M — fl, by calculating 
the elliptic modulus of the semiclassical solution. 

For g — 0, Eq. I|A7(1 describe all possible dynamics 
for arbitrary detuning energy u), while in the presence 
of a pairing interaction the dynamics is given by solu- 
tions l|A15(l and l|A19|) . depending on the sign of D. The 
dynamics is described by Eq. (|A20I) at a single singular 
point in k — r\ space, as we shall see later. 

Figure ITU shows the elliptic modulus of the semiclassi- 
cal solution in k — rj space. The two regions D > and 
D < are separated by D = lines that correspond to 
the specific values of elliptic modulus k = and k = 1. 
Eqs. I|A15|) and (|A19|) coincide for k = 0, which corre- 
sponds to the white lines in Fig. Hence these two 
solutions connect continuously when crossing that line. 
For the black line k = 1, on the other hand, solutions 



FIG. 14: Elliptic modulus in k — g space for the half-filled 
shell: 2M = 



discontinuity gives rise to the self-trapping transition dis- 
cussed in subsection IIV CI 

Figure [21 shows that the black line and one of the 
white lines intersect at the critical point (k c ,i) c ), where 
q = p = 0, given explicitly by 



Vc 



2^1(45 + 26^3) 1/4 ~ 5.0302 (A21) 

(A22) 



k c = (2 + V3) ~ -2.9677 



Form this result, we conclude that the self-trapping tran- 
sition appears by varying the detuning parameter k only 
for 77 > rj c . 
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